An Itzykson-Zuber-like Integral and Diffusion for Complex Ordinary and 

Supermatrices 



> 
o 



in 



Thomas Guhr a > and Tilo Wettig b > 
Max Planck Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 

(May 15, 1996) 

We compute an analogue of the Itzykson-Zuber integral for the case of arbitrary complex matrices. 
The calculation is done for both ordinary and supermatrices by transferring the Itzykson-Zuber 
diffusion equation method to the space of arbitrary complex matrices. The integral is of interest for 
applications in Quantum Chromodynamics and the theory of two-dimensional Quantum Gravity. 

PACS numbers: 02.30.Cj, 11.30.Pb 



I. INTRODUCTION 



In 1980, Itzykson and ZubeiE presented their result on a certain integral over the unitary group which had great 
impact in several areas of mathematical physics. Let U be a matrix parameterizing the unitary group U(N) with the 
invariant Haar measure dfi(U). Moreover, consider two diagonal matrices x and y with entries x n and y n , respectively, 
I/"") i where n — 1, . . . , N. The Itzykson-Zuber integral can then be written in the form 



dn{U) exp(ztrf7 l xUy) = , \ > L1 

A N (x)A N (y) 



where the function 



^N(x) = Y[ {Xn-X m ) (1.2) 

is the Vandermonde determinant of order AL Although it was later realized that this formula is a special case of a 
more general result due to Harish-Chandra,El it prompted many investigations in various fields. In 1983, Mehta and 
Pandejola used this formula to work out, in the framework of Random Matrix Theory, the spectral correlations of 
a generic quantum chaotic system which undergoes a transition from conserved to broken time-reversal invariance. 
There are also numerous applications in field theory, particularly in the theory of two-dimensional Quantum Gravity; 
^ (— | a review can be found in Ref. |[ 

Recently, ShatashvilitJ showed that the integral (1.1) itself and all correlations in the Itzykson-Zuber model can be 



evaluated using the Gelfand-Tzetlin coordinates for an explicit calculation. Remarkably, Itzykson and Zuber had not 
derived their result by an explicit calculation but related the integral ([0]) to a diffusion process. They showed that 
it can be viewed as the kernel of a diffusion equation in the curved space of the eigenvalues of Hermitean matrices. 
Since the space of Hermitean matrices is Cartesian and, therefore, easy to treat, the result ([Pi]) can be found by 
comparison with the curved space without actually calculating it explicitly. The crucial point is the separability of 
the Laplacian operator in the curved space of the eigenvalues. This diffusion equation technique turned out to be a 
very powerful tool. _ 

Some years ago, it was realized by one of the present authorso that the Itzykson-Zuber diffusion can be directly 
transferred to supermathematics. After the pioneering, mathematical achievements of Berezinjj superrnathematics 
was brought into the theory of matrix models by EfetovO and Verbaarschot, Weidenmuller, and Zirnbauerli!] In Ref. 
it was shown that there is, completely analogously to the case of ordinary matrices, a diffusion equation in the space of 
the eigenvalues of Hermitean supermatrices whose kernel is the supersymmetric generalization of the Itzykson-Zuber 
integral (1.1). Thus, the integral can be worked out generalizing the methods of Ref. |. Again, the crucial point is 



the separability of the Laplacian in the curved space of the eigenvalues E The result has been used to evaluate, in 
the framework of Random Matrix Theory, the effect of symmetry breaking on the spectral correlations of a chaotic 
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time-reversal non-invariant system.0 Recently, it was shown that the supersymmetric generalization of the Itzykson- 
Zuber diffusion has not only a mathematical, but also, a physical meaning: in Random Matrix Theory, it describes 
the transition from arbitrary to chaotic fluctuations .113 These results could be used to work out the crossover from 
Poisson regularity to chaos in a time-reversal non-invariant system. 

Let u be a supermatrix, parameterizing the unitary supergroup U(ki/k2), with invariant Haar-Berezin measure 
and let s and r be diagonal matrices, both having the form s = diag(si, is^ with Sj — diag(sj r j ir _ r . , Sfcj), 
j = 1, 2. Then, in its most general form, the supersymmetric Itzykson-Zuber integralQ can be written asl 13 tl 14 l 

M u) ex P (itrg^-V) = - ^<"^> W=i.-... fcl p i^p^ksl^ , (1 . 3 ) 

B kl k 2 {s)B klk2 {r) 

where the symbol trg stands for the supertrace, often also denoted by str. The supersymmetric generalization of the 
Vandermonde determinant is given by 

B klk2 {s) = Ti : -. (1.4) 

Up,g{ S pl - lS q2) 

It is worth mentioning that B kl k 2 ( s ) reduces to a determinant^ via Cauchy's lemma in the case fci = k 2 such 

that -Bfc, fcj (s) = det[l/(s p i — is a 2)] v , a =x k, ■ Furthermore, by setting kx = N and k 2 = 0, one easily sees that 

formula ( |l.3| ) includes formula (1.1) as desired. However, for k\ and k 2 both non-zero, there is an important caveat: 
In superanalysis, a change of variables can ind uce a certain class of singularities. Here this implies that, if further 
integration over, say, the s- variables in Eq. ( |l.3[ ) is required, we have to deal with new types of boundary contributions 
which have no analogue in ordinary-, analysis. The general theory of this effect, which is sometimes overlooked in the 
literature, was given by Rothstein.Ej In Refs. and [l6| it was shown how to treat these contributions in the case of the 
supersymmetric Itzykson-Zuber integral. Regarding integrations over-. the. sup ersymmetric Itzykson-Zuber integral , ye t 



ano the r comment is in order: in applications in statistical mechanics, BuLrEilO the integrand, containing formulae (1.1) 



EO C3 



or (1.3), possesses some invariances under permutations. This allows one to replace the numerator of the right hand 
sides of these equations by exp(itvxy) or exp(itrgsr), respectively, which makes the calculations more transparent. 
We emphasize this trivial point since it has stirred some confusion lately. 

It is the purpose of this work to transfer all the results which have been reviewed so far from Hermitean to 
arbitrary complex matrices. We will do this for the case of ordinary and supermatrices. We will derive closed 
formulae for the analogues of the ordinary and the supersymmetric Itzykson-Zuber integral in the space of complex 
matrices. To the best of our knowledge, those have not been worked out yet. Our results are of considerable interest 
in the theory of ma t rbt .models. Several models in two-dimensional Quantum Gravity involve complex instead of 
Hermitean matrices. LZJL3 Recently, the so-called chiral Gaussian Ensembles,t2l also based on complex matrices, 
been introduced and proved to be very useful in studying certain aspects of Quantum Chromodynamics (QCD).I 
In particular, the integral we compute will be crucial for a further analysis of the spectral correlations of the Dirac 
operator of QCD in the framework of random-matrix models atpfinite temperatures, in the presence of mass terms, 
or at finite chemical potential. Via the Banks-Casher formula,L3 such analysis is also important for the study of 
questions related to the chiral phase transition of QCD. 

To derive our results, we transfer the diffusion equation technique to complex ordinary and supermatrices. Again, 
the crucial point turns out to be the separability in the curved space of, in this case, radial coordinates. We have 
no doubts that the explicit formula could also be evaluated using other techniques. Especially, Gelfand-Tzetlin 
coordinates could be used as in Ref.-g] for ordinary Hermitean matrices. Recently, Gelfand-Tzetlin coordinates were 
derived for Hermitean supermatriceaU such that this method could also be used for complex supermatrices. Moreover, 
regarding the case of ordinary matrices, we do not exclude the possibility that our results might be derivable directly 
from Harish-Chandra's formula.B To the best of our knowledge, Harish-Chandra's result has not been transferred to 
superanalysis yet. Nevertheless, here we will focus on the diffusion equation technique. So far, this technique was, for 
Hermitean matrices, viewed as a purely mathematical tool. However, as mentioned beforftr-the diffusion in superspace 
describes the transition from arbitrary to chaotic fluctuations in Random Matrix Theorylij and, therefore, also has a 
direct physical meaning. We strongly believe that similar features are likely to exist in the case of the diffusion in the 
space of complex matrices which we will discuss in the present work. 

The paper is organized as follows. In Sec. |J, we state our results and derive them by constructing an eigenvalue 



equation. In Sec. Ill, we discuss the diffusion and questions related to it. We summarize and discuss our results 



in Sec. [V. Three appendices are provided for the derivation of intermediate results used in the text. In a fourth 



appendix, we discuss some boundary contributions which occur in the case of Hermitean Supermatrices. 
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II. DERIVATION OF THE INTEGRAL BY CONSTRUCTING AN EIGENVALUE EQUATION 



In Sec. II A, we state the integral for both case s, for ordinary and supermatrices. The derivation is performed for 
ordinary and supermatrices in Sees. II B and II C, respectively. 



A. Statement of the Integral 

Let X be an arbitrary square, complex ordinary matrix of dimension N. It is well knownB that it can be written 
in the so-called pseudo-diagonal form 

X = UxV with x — diag(a;i, . . . , xn) , (2.1) 

where the N variables x n are real and non-negative. Note that these are not eigenvalues, they will be referred to as 
radial coordinates. While the matrix U explores the full parameter space of the unitary group U(iV), the matrix V is, 
in order to remove a double counting of phases, restricted to a subspace defined as the quotient of the unitary group 
and the Cartan subgroup, hence we have U £ U(JV) and V £ U(iV)/U JV (l). 

We now multiply A" by a diagonal matrix y of the same form as x and consider the expression iRetvUxVy. The 
integral we wish to compute is the angular average over both unitary matrices of the exponential of this trace, 

$(x,y) = J dfj,(U) J dfi(V) exp(iRetrUxVy) , (2.2) 

with dfi(U) and d/i(V) being the corresponding invariant Haar measures. We show that this integral is given by 

*/ v (27r) w2 det[J (x„y ro )] x ;v 

* (:C ' y) — A^A^) ' (2 - 3) 

where Jq(z) is the ordinary Bessel function of zeroth order. The Vandermonde determinant was defined in Eq. ( [1.2] ). 

Remarkably and fully analogously to the case of Hermitean matrices, this result can straightforwardly be generalized 
to supermatrices. An arbitrary complex supermatrix a of dimension ki + &2 can be written as 

a = usv with s = diag(si, isz) (2-4) 

and Sj — diag(sij, . . . , Sk 6 j) for j = 1, 2. Again, while the matrix u explores the full parameter space of the unitary 
supergroup, the matrix v has to be restricted to a subspace in order to remove phases, we thus have u £ U(fci/A;2) 
and?; G U(fei/fc 2 )/U i:i + fc2 (l). 

Analogously to the case of ordinary matrices, we multiply the matrix <r by a diagonal matrix r of the same form as 
s and consider the expression Re trgusvr. The generalization of the integral in the ordinary case is given by replacing 
the trace by the supertrace and the invariant measures by the corresponding ones in superspace, dfi(u) and dfi(v), 
respectively. We show that the double average 



f{s,r) = J dfj,(u) J dji(v) cxp (i Re trg usvr) (2-5) 

is given by 

_ (2^)( fc i~ fc 2) 2 det [Jo(a P iyi)] p , p , =1 ,... [fcl det [Jota^r^)]^ fa 

nS ' r) 22*fafc 1 !fc 2 ! B klk2 (s 2 )B klk2 (r 2 ) ' 1 ' 



where the generalized Van dermonde determi nan t was defined in Eq. (1.4). 



Obviously formula (2.3) includes formula (2.6) as can be seen by putting k\ = N and &2 = 0. Hence, in principle, it 
is sufficient to perform the derivation solely in superspace. However, we decided not to do so. We prove both results 
separately, first, in order to give those readers with little interest in supermathematics the opportunity to understand 
the ordinary case without being burdened by undesired information and, second, in order to help those readers with 
little experience in supermathematics to approach this area starting from more familiar grounds. 

Note that our discussion is related to the harmonic analysis in the corresponding matrix spaces. The functions 
3>(x,y) and ip(s,r) can be viewed as the lowest order Bessel functions in these spaces. 
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B. Derivation for Ordinary Matrices 



Besides X, we introduce a second arb itrary complex matrix Y of dimension N whose pseudo-diagonalization reads 
Y = U'yV with y defined in Sec. |LA|. We have U' £ U(JV) and V' £ U(N)/U N (1). We observe that the "plane 
wave" 

W{X, Y) = exp(iRetrXy f ) (2.7) 
in this matrix space obeys the "wave equation" 

AW(X, Y) = — (trFF f ) W(X, Y) , (2.8) 

where the Laplacian is defined as 

v / d 2 d 2 \ 



u rn ) 



Due to the invariance of the Haar measures we can express the function ( |2.2j ) as the angular average of the matrix 
plane wave, 

${x,y) - / dfi(U') ( dfi(V')W(X,Y) . (2.10) 



The crucial observation is now, just as in the case of Hermitean matrices, that &(x, y) satisfies the wave equation 
This can be seen by averaging both sides of (2.S) over the matrices U' and V' using trFyt — try 2 . Consequently, since 
4>(x,y) depends only on the radial coordinates, we can replace the Laplacian A by its radial part A x . To construct 
it we have to transform the Cartesean volume element 



to radial and angular coordinates, 



d[X] = Y[ dReX nm dImX nm (2.11) 



d[X] = J(x)d[x}dn(U)d/j(V) 

N 



d[x] = J^J dx r , 

n=l 

N 

J(x) = A 2 N (x 2 )Y[x n , (2.12) 



where the Jacobian J(x) was worked out in Ref. It is then easily shown that the radial part of the Laplace operator 
reads 

N 1 d d 

^-j J(x) dx n dx n ' ^ ^ 

n— 1 x ' 

and we thus arrive at the eigenvalue equation 

A x ${x,y) = -{try 2 )<S>{x,y) (2.14) 

in the curved space of the radial coordinates. 

The key to the solution of the above equation is the separability of the radial Laplacian. For an arbitrary function 
A(x) we have the identity 

A N (x 2 ) A N (x 2 ) 
where the reduced part of the radial Laplacian is 
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A' 



A * = E 



1 d 

\ dx^ x n Qxji 

The proof of this fact is given in Appendix |X| Hence, the ansatz 



A w (x 2 )A JV (y2) 



(2.16) 



(2.17) 



in which, for symmetry reasons, x and y are treated on the same footing, reduces the radial equation (2.14) to the 
much simpler form 



A' x y(x,y) = ~{try 2 )y(x,y) 



(2.18) 



This equation is again separable by a product ansatz for ^(x, y) which yields TV Bessel differential equations of zeroth 
order for each of the N functions. Hence, each of them can be written as a linear combination of the Bessel and Weber 
functions Jo and A*o, respectively. Note that the indices of the eigenvalues y\ of these N Bessel differential equations 
can be permuted arbitrarily implying that the most general solution is a linear combination of all these permuted 
products. However, the integral representation (2.2) imposes certain boundary conditions on the solution of the 
differential equation ( |2.1g| ). Since the integral has a finite value for all x and y we have to exclude the Weber function 
from the solution. Additionally, we have to take into account that the integral is invariant under permutations of the 
indices. Since An(x 2 ), the Vandermonde determinant, is really a determinant, the function ^>(x,y) has to have the 
same properties under permutations of the indices. Incorporating these boundary conditions we find 



(2tt) n2 

®(x,y) = , det [J (x n y m )] r 



,N ' 



(2.19) 



which yields immediately the result (2J3). Of course, the normalization constant is arbitrary. We will show later why 
our choice is useful. 



C. Derivation for Supermatrices 



All steps are completely analogous to the ordinaW|-&ase. In order to make the notation more transparent, we write 
the supermatrix a in the boson-fermion block forrroliJ 



(2.20) 



where cr J ' is a kj x ki complex matrix whose entries are commuting if j = I and anticommuting if j ^ I. The factor 
i in front of a 22 is, as usual, introduced to ensure convergence. ErE3 Again, besides a, we introduce a second arbitrary 
complex supermatrix p of the same form whose pseudo-diagonalization reads p — u'rv 1 with r defined in Sec. [I A and 
with u' € U(Jfei/& 2 ) and v' e V(k 1 /k 2 )/\J kl+k2 (l). There is also a "plane wave" 



w(a,p) = exp(zRe trg ap^) 



(2.21) 



in this matrix space. Note that the expression Retrgcrp^ has, for explicit calculations, to be interpreted as half the 
sum of trgcrpt and its complex conjugate since we will not introduce the real and the imaginary part of Grassmann 
variables. The plane wave satisfies the "wave equation" 



Aw(a, p)=- (trgpp 1 ) w(a, p) 



where the Laplacian is defined as 

2 



t) 2 



4 EE 

3^1 P. 1 



BfT jl rr 3 '* 



(2.22) 



(2.23) 



The invariance of the Haar measures allows us to express the function (2.5) as the angular average of the matrix plane 
wave, 
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<p{a,r) = / dfi(u') / dn(v')w(a,p) . (2.24) 



As in the ordinary case, we integrate both sides of t he w ave equation over the matrices v! and v' using trgppt = trgr 2 
and observe that (p(s, r) satisfies the wave equation ( 2.22| ). Again, since (p(s, r) depends only on the radial coordinates, 
we can replace the Laplacian A by its radial part A s . The transformation of the Cartesean volume element 

2 

d[*] - n n d ( Re < 9 ) rf ( im <j n n ^ 

j'=i p,q jjti P<<i 

to radial and angular coordinates reads 

d[cr] = J(s)d[s]dfjt(u)dfj,(v) 

2 kj 

d ^ = n n ds Pi 
j=i P =i 

2 kj 

J(s)=B 2 klk2 (s 2 )l[ l[s pj , (2.26) 

j=i p=i 

where the Jacobian or Berezinian J{s) is computed in Appendix |^. The radial part of the Laplace operator takes the 
form 



2 k 

details are given in Appendix |b| Hence, we have to solve the equation 

A s <p(s,r) = -(trgr 2 )(p(s,r) (2.28) 

in the curved space of the radial coordinates. 

In the case of ordinary matrices, the key for the solution was the separability of the radial Laplacian. It is essential 
that this_feature is also present in the case of supermatrices. This closely parallels the situation for Hermitean 
matricesfl For an arbitrary function A(s) we find 

s AM ' A^A( S ), (2.29) 



where the reduced part of the Laplacian reads 

2 kj 



ai-EElBT+rrsH ■ (2 - 30) 

j=lp=l \ OS P3 S P3 OS P3 ) 



The derivation is given in Appendix ^|. Thus, the ansatz 

B klka (s*)Lk l i,Ar-i 



yields the reduced equation 

A'>(s,r) = ~(trgr 2 )ifj(s,r) (2.32) 

which, again, is separable by a product ansatz for ip(s,r). We obtain k\ + It2 Bessel differential equations of zeroth 
order. The boundary conditions imposed by the integral representation (2.5) are very similar to the ones in the 
ordinary case. First, we have to construct the solution using the Bessel function Jo and to reject the Weber function 
Nq. Second, we have to take into account the invariance under permutations. Here, however, we see from the Jacobian 
that this invariance exists only within the boson-boson or fermion-fermion block, respectively. Since these boundary 
conditions imply that the solution is given by 

(27r)( fe i- fe 2) 2 

^ r) = 2 2fc Afci!fc 2 ! det [ J °( s P ir P'i)] P ^=i,...,fc 1 det [^o(s,2V2)]^ = i fe , (2.33) 

we arrive at the final result (2.6). Again, the normalization constant is arbitrary, we will comment on our choice later. 
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III. DIFFUSION EQUATION AND FOURIER TRANSFORM 



We now discuss a diffusion equation which is closely related to the plane waves and the eigenvalue equations we 
constructed in the previous section. Our goal is to show that the concept of diffusion equations which is so useful in 
the case of Hermitean matrices can be transferred straightforwardly to arbitrary complex matrices. However, since 
these considerations are more of conceptual interest and do not require so many explicit calculations, we study only 
the case o f supermatrices. The case of ordinary matrices is always recovered by setting k\ = N and ki = 0. In 



Sec. niB 



Sec. |III A , we intr oduce the concepts in Cartesian space. We go over to the curved space of the radial coordinates in 



In Sec. [II C, we discuss some questions related to boundary terms. 



A. Cartesian Space 



We introduce a time coordinate t and consider the diffusion equation 



for a given initial condition F (a) such that 



\imF(a,t) = F (a) . 



The kernel of this diffusion satisfies the equations 



\ AG(<7, t) = ^-G(a, t) and lim G(a, t) = 5(a) , 



(3.1) 



(3.2) 



(3.3) 



where the 5-function is given by 



w = n n <) s ( im <) n n 



J ' ) 



(3.4) 



The 5-function of an anticommuting variable (3 is defined by 5(/3) = VSTr/^-Q'^l Similar to the discussion in Refs. |l| 
and |l^, the kernel is the Gaussian 



G(a,t) 



exp ( -— trgaa 1 



(27Tt)( fc l-^) 2 

and the solution of the diffusion process can be written as the convolution 

F{a,t) = J G(a-a',t)F (a')d[a') . 
Moreover, this solution is also expressible as 



F(a,t) = exp ( -A ) F (a) 



(3.5) 



(3.6) 



(3.7) 



which has to be viewed as a formal power series. 

We will now show how the diffusion can be related to the pla ne w aves of the previous section and to the theo ry of 
Fourier transforms. To this end, we remark that the 5-function (3.4) can be expanded in the plane waves (2.21), 



5(a) 



w(a, p)d[p] 



( 27r )2(fc 1 -fc 2 )2 

which allows us to introduce the Fourier transform P(p) of a function P(a) and its inverse by 



P(P) 



P(a)w*(a,p)d[a] and P(a) 



22/cife 2 
( 27r )(fcl-fc 2 ) 2 



P(p)w(a,p)d[p] . 



(3.8) 



(3.9) 
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The Fourier transform can be used to derive the explicit form (3.5) of the diffusion kernel defined in Eq. (3.3), this 
works as follows. The diffusion kernel can, according to Eq. (3/7), be expressed in the form 



G(*,t)=exp(-A)8{a) 



(3.10) 



in which we insert the expansion (3. 



G(a,t) 



24fcifc 2 



(2 7r )2(fe 1 -fc 2 ) 2 



exp ( - A ) w{a,p)d[p] 



(3.11) 



We write the exponential as a power series and, by virtue of the eigenvalue equation ( |2.22 ), perform all derivatives. 
The resummation of the series gives the diffusion kernel as the Fourier transform of a Gaussian 



G{<7,t) 



( 27r )2(fci-fc 2 )2 



exp 



trgpp f ) w(a,p)d[p] , 



(3.12) 



which is in agreement with Eq. (3.5). 



B. Curved Space of the Radial Coordinates 



We now assume that the in itial condition depends only on the rad ial coordinates, i.e. Fq{ct) — Fq(s). Thus, it is 
useful to use the coordinates (2.4) in the integral representation ( |3.6[ ) of the solution of the diffusion equation (|3.l|). 
This has some important consequences. The invariance of the Haar measures implies that this solution is also a 
function of the radial coordinates only, hence we have F(a, t) — F(s, t). Consequently, the diffusion takes place in the 
curved space of the radial coordinates alone, 



X - A S F( S , t) = t) and Urn F(s, t) = F (s) , 



(3.13) 



where A s is th e ra dial part of the Laplacian define d in Eq. (2.27). Moreover, we may conclude from the integral 
representation (3J3) that the kernel of the diffusion ( |3.13] ) is given by 

T(s,s',t) = J dp{u) J dp(v)G{usv- s',t) . (3.14) 

There are t wo ways of e valuating this double average. First, since the kernel G(a,t) is Gaussian, a direct comparison 
of Eq. (|J) to Eq. (|J) shows that 



T(s, S ',t) 



( 27rt )(fcl-fc2) 2 



exp 



(-ltrg( S 2 + ,s' 2 )) V (-isM 



(3.15) 



which means that this double averag e ca n be expressed in terms of the one we have calculated in the previous section. 
Hence, with the help of the result (2.6) and after a reordering of factors, we can write the diffusion kernel in the 
curved space in the form 



r(s,s',t) 



1 d(3t b( S Pl' S p'l^)] p ^ = 1: ... !fcl det [7(« 9 2. S 'q'2^)] q>q , = h 



ki\k 2 \ B klk2 ( S 2)B klk2 (s>2) 
in which the entries of the determinants are given by the function 

1 



j(s pj ,s' qj ,t) = - exp 



2 1 12 



21 



t 



(3.16) 



(3.17) 



for all values of j — 1, 2 and p, q — 1, .. . , kj. The function Ip(z ) is the modified Bessel function of zeroth order. 

Alternatively, if the result ( |2.6| ) was unknown, formula (3.16) could be derived by a procedure similar to the one in 
Sec. IIC. The separ ability of the radial part A s leads to a reduced diffusion equation involving the reduced operator 



A' s defined in (2.30). This equation can be solved by a product ansatz leading to the diffusion equation 



8 



1 / d 2 



2 \ d: 



PJ 



1 d_ 



d 

l(s pj ,s qj ,t) = -Q t l(s PJ ,s' qjl t) 



(3.18) 



where the differential operator is just the radial part of the Laplacian in a two-dimensional space. In order to construct 
the solution, we express it as the formal series 



l{s pj ,s' ,,t) = exp - 



t / d 2 



1 d 



S PJ S qj 



(3.19) 



acting on the proper radial (5-function in this two-dimensional space. Inserting Hankel's expansionEl of this ^-function, 

J (Spjz)J {s' qj z)zdz , (3.20) 



3( S Pj S qj, 



we can perform all derivatives and resum the series. We arrive at 



oo 

l{ S PV S 'qj^) = / eXP (~2" Z2 ) J °( S P3 Z ) J o( S 'qj Z ) zdz ' 



(3.21) 



which is precisely Weber's representationEj of the function (3.17). It is easy to see in a direct calculation that this 
function is indeed the kernel of the diffusion equation ( 3.18 ). According to an elementary resulto of the theory of 
Bessel functions, 7(s PJ , s ' q j,t) is properly normalized, 



oo 

/• 



"f(s p j, S q j, t) SpjdSpj — 1 



(3.22) 



where we have used the radial part s P jds P j of the measure in the two-dimensional space. Furthermore, since Io(z) 
behaves like exp(z)/v / 2~7rz for large values of the argument, the function "/(s P j,r q j,t) approaches the 5-function 



lim 7 {s pj ,s' ,,t) = 



5(s p j s„j) 



(3.23) 



S PJ S qj 



for vanishing time t. 

The limit relation (3.23) implies that the kernel ( |3.16 ) satisfies the correct limit relation in the curved space of all 
radial coordinates, we write this in the form 



lim r(s, s' , t) = I dfi(u) / d/i(v) 5(usv — s') 



klW- y/J(s)J( S ') 



(3.24) 



where the Berezinian J(s) is defined in Eq. ( 2.26 ). Using this result, it is easily checked that the constant l/fcjjfcg! 
ensu res th e correct normalization. This, in turn, motivates our choice of the normalization constants in Eqs. ( 2.19| ) 
and (|2~33|). 



C. Questions Related to Boundary Contributions 



The function T(s,s',t) given in Eq. (3.16) is, as we have shown, the kern el of the diffusion equation (3.13) in 



the curve d sp ace of the radial coordinates. Thus, the solution of the integral (3.14) as it stands is indeed given by 
formula ( 3. 16 ). However, there is a very subtle point about kernels of this type which has an important impact 
on applications. Although we will present some applications of our results to physical problems in a forthcoming 



9 



publication, we already give a short, more intuitive, discussion of this subtlety here. We do so to acquaint the reader 
who is not yet familiar with supersymmetry with this point which is som etim es overlooked in the literature. 
In Cartesian space, the normalization of the Gaussian diffusion kernel (3.5) implies that the equation 



G(a-a',t)d[a'\ = 1 



(3.25) 



holds for all values of t and for all matrices a. Thus, after performing the angular integration, one would naively 
assume that the radial integral 

rj(s,t)= Jr(s,s',t)J(s')d[s'] (3.26) 

also yields unity for all values of t and for all diagonal matrices s. In the ordinary case, i.e. for k\ = N and k2 = 0, 
it can be checked by a straightforward, direct calculation that we indeed have r](s,t) = 1. However, in the case of 
supermatrices this is, for non-trivial reasons, no longer true. The singularities of the Berezinian J(s) compensate the 
vanishing of some angular Grassmann integrals such that a finite, non-zero result remains. This effect leads to certain 
contributions to the integral which are often called Efetov-Wegner-Parisi-Sourlas terms in the more physics-Qdexi±ed 
literature. There are various methods to construct those contributions in the applications of supersymmetryElalij'EZl 
From a strictly mathematical point of view, these terms arise as boundary contributions due to the fact that the 
integrals are ill-defined for non-compact supermanifolds. A full-fledged theory can be found in Ref. [l5L It is instructive 
to think of these boundary contrib ution s as being necessary to r estor e the translational invariancp,of the integrals in 
Cartesian space, as obvious in Eq. ( 3.25| ), which is broken in Eq. (3.26) if r/(s, t) differs from unity.c3 To illustrate this, 
we calculate the function ry(s,t) for the simplest non-trivial case, namely k\ = k 2 = 1. We have 



r)(s,t) = (s 



OO CO 

2 4- « 2 \ f f di' ds' S ' llS ' 12 
n + s i2> J J as n as 12 J ^ ^ 



j(sil,s' n ,t)j(si2,s' 12 ,t) 



(3.27) 



By expressing the denominator of the Berezinian as the integral 

OO 

1 



exp (-(si 2 + s' 2 2 )z) dz 



(3.28) 



we can evaluate the double integral (3.27) by standard methodsEj We arrive at 

r](s, t) = 1 — exp 



s ll + s 12 



2t 



(3.29) 



which equals unity only in the limit t — > 0. Note that the exponential is, apart from a numerical factor, nothing else 
but the Cartesian kernel at a' — which is just G(a, t) = G(s, t). This is, of course, no accidental coincidence. 

The case fci = where k\ is arbitrary, is physically the most interesting one. Due to the deter minant structure of 
the function Bk 1 k 1 (s' 2 ), the evaluation of the function r](s,t) reduces to the double integral ( 3.27 ) already computed, 
and we arrive at 



rj(s,t) 



1 



B klkl (s 2 ) 



det 



1 — exp 



'pi 



S pl + S q2 
2t 



i.q— l,...,ki 



As evident from the definition (3.26), this function is a solution of diffusion equation (3.13), we have 



1 . . d 
-A g r){s,t) = qIVKS,*) ■ 

However, it is not a kernel of the diffusion process in the usual sense since it obeys different limit relations, 



lim n(s, t) 
t-»o 



1 and lim rj(s, i) = , 



(3.30) 



(3.31) 



(3.32) 



which reflect the existence of the new boundary contributions. The function rj(s, t) can be viewed as the envelope 
solution corresponding to the kernel T(s,s',t). The function Bk 1 k 1 {s 2 )T(s, s' ,t) possesses a product structure in the 



10 



kernels 7(s p i, s' 2 , t). Due to the integration, this property has vanished in Bk 1 k 1 (s 2 )r](s,t), which factorizes only in 
functions of the combinations + s^ 2 - Along the lines given in Refs. [j], [l^, and |2^, one can show that the diffusion 
kernel in the curved space of the radial coordinates has to be replaced by 

T(s, at, t) — (1 - V(s, t))^r + T(s, s\ t) (3.33) 

if further integration over the primed variables s' is required. This replacement cures the problem of the boundary 
contributions for k\ = k 2 in the physically most interesting cases. Hence, for a well-behaved initial condition Fo(s), 
the solution F (s, t) of the diffusion equation in the curved space of the radial coordinates reads 

F(s,t) = (l-77(a,t))F (0) + / ' T(s, s' ,t) F (s') J(s')d[s'} . (3.34) 

We emphasize that this result is really a solution of the diffusion process ( |3.13| ), including its initial condition. Note 
that there are very peculiar cases in which further boundary contributions can arise. Those, however, have to be 
constructed using the full theory which is given in Ref. ^|. 

In Appendix [pj, we reconsider the boundary contributions to the supersymmetric Itzykson-Zuber integral for Hcr- 
mitean Supermatrices as derived in Ref. [7l 



IV. SUMMARY AND DISCUSSION 



We have calculated an analogue of the Itzykson-Zuber integral in the space of arbitrary complex matrices. We 
arrived at explicit formulae for the case of ordinary and supermatrices, where the latter includes the former. We 
performed our calculation by transferring the diffusiatL-equation technique of Itzykson and Zuber for Hermitian 
matrices, which works in ordinarju and in superspacejlltj to complex matrices. For the actual derivation, we used 
an eigenvalue equation for the plane waves in these matrix spaces which is closely related to this diffusion. Similar 
to the Hermitean case, the integral in question turns out to be the kernel of the diffusion in the curved spaces of the 
radial coordinates. The explicit results can be computed due to a separability of the Laplacian in these radial spaces. 
We discussed certain types of boundary contributions to the full solution of the diffusion equation which can arise in 
superspace. 

We have no doubts that our explicit results can alsa,be derived by other methods. In particular, the use of Gelfand- 
Tzetlin coordinates for the unitary group in ordinarya and superspaceE-3 ought to be mentioned here since it allows a 
recursive evaluation of correlation functions in the corresponding matrix models. Most importantly, as far as the case 
of ordinary matrices is concerned, it does not seem unlikely that the explicit formula for the integral we presented 
here can be derived directly from the Harish- Chandra integral. At first sight, one would not think so since the 
Harish-Chandra integral is an average over one unitary matrix whereas our result is a double average over two unitary 
matrices. This can be seen from the fact that our explicit formula contains Bessel functions where the Itzykson-Zuber 
integrals contain exponentials, i.e. plane waves. The Bessel function of zeroth order is just the angular average over 
a plane wave in a two-dimensional space. This might indicate that the angular average over two unitary matrices 
is essential and can not be replaced by an average over one unitary matrix. Nevertheless, we do not exclude the 
possibility that a clever reordering of the trace in the matrix plane waves which combines these two unitary matrices 
can be done in such a way that the essential part of the calculation reduces to an application of the Harish-Chandra 
formula. These considerations, however, do not apply to the case of supermatrices, for which, to the best of our 
knowledge, Harish-Chandra's result has not been transferred yet. 

For the reasons discussed above, we do not want to present our explicit formulae for the integrals as our most 
important findings. Rather, we consider the Itzykson-Zuber-like diffusion which we constructed here as our most 
interesting result. We strongly believe that this diffusion is more than a mathematical tool to calculate integrals. In 
the case of Hermitean supermatrices, it was shownlij that the Itzykson-Zuber diffusion models the transition from 
arbitrary to chaotic fluctuations of all orders in a very general way. We are convinced that the diffusion in the space 
of complex matrices also has a physical meaning of similar significance. 
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APPENDIX A: SEPARABILITY IN ORDINARY SPACE 



We use a more convenient form of A x , 



Ax - J2[ 



d 2 d\nJ(x) d 



dxi dx n dx n 



(Al) 



The derivatives in Eq. ( [2.15 ) are evaluated in a straightforward manner, and we arrive at the intermediate result 



Ax 



1 N 



d 2 1 d 



A N (x 2 ) A N (x 2 ) ^-J \dx 2 t ' x n dx n 



- 4D n A(x 



(A2) 



where D n = S n + x 2 n (S 2 - T„) with 



N 1 



and T n - ^ - . 

m= i k^n ^m) 



(A3) 



We now show that Xm=i ^« = 0- The ^ ac * ^ na ^ 



n— 1 m^tn '* "* 



(A4) 



is easily seen by renaming summation indices. The remaining term is 

2V 



/'/y>2 . y 2 *\ / rv* 2 ____ \ 



n=l 



S2( 



(A5) 



where the symbol f2(n, m, m') denotes summation over three indices n,m,m' which are pairwise different. We now 

—R , (A6) 



2R= £ 



(rv»2 y» 2 ^ z' /->i 2 o^2 \ / /r»2 *\ /' q-* 2 -vi 2 ^ 

x m J 'n/V' i 'm x mV V x m' • L m)\ J 'm' • L n) 



from which i? = follows immediately. This completes the proof. 



APPENDIX B: DERIVATION OF THE BEREZINIAN AND THE RADIAL PART OF THE LAPLACIAN 



We wish to compute the Berezinian of the transformation a — usv, defined in Eq. (2.4). We first construct 

da = u(w dus + ds + sdvv^)v . (Bl) 

Writing v)du = du' and dvv* = dv' , and noting that s^ — s, we obtain the invariant length element 

trg dada^ — trg (du s + ds + sdv) (sdu 1 * + ds + efii'^s) (B2) 

from which the Berezinian can be read off. Since du' and dv' are also in the algebra we are entitled to drop the primes. 
This gives 



trg dada^ = trg ds 2 + trg sds (du + dv) + dv + e?i>^) + trg (dus + sdv) (sdu' + dv^ s) , 



(B3) 



where we have made use of the fact that s and ds are diagonal. Since du and dv are anti-hermitian, the second 
term in the above expression yields zero. The vanishing of this term also shows that the Laplace ope rator separates 
into two sums over radial and angular coordinates, respectively, a fact which has been used in Sec. II C. The first 
term in Eq. ( ]B3| ) contributes a factor of one to the Berezinian so that we are left with the third term only. Using 
boson-fcrmion block notation, we write 



12 



du 



du c 



-du At du° 2 



(B4) 



A note about the number of degrees of freedom: du j is an anti-hermitian matrix with k 2 commuting degrees of 

freedom whereas du A and du A ^ each have k\k 2 anticommuting degrees of freedom. Similar notation is used for dv, the 
main difference being that the diagonal elements of dv Cj are zero. It is convenient to separate the diagonal elements 
of du Cj and to define 



dus + sdv = rj + u 



' y 1 







' UJ 11 UJ 12 ' 





if 2 


+ 


UJ 21 UJ 22 



where rj 11 = diag(duf 1 1 sn, 
u! 22 are zero, and 



'dUhk^kii), rf 2 = diag(dwf 1 2 isi 



(B5) 



du^ k<i isk 2 2), the diagonal elements of uj 11 and 



i dv. 



Ci 



m 

UJ 



u> 



pi 

21 

gp 



qq' 

12 - du A q is q2 + s pl dv A q 



pp 



ip + V) 



du qp Spi - is„idv„„ 



We now consider 



trg (rj + uj) (?7 t + cjt) = trg (?777 t + ujuj^) 



= trg W t + tr uj 11 ^ - tr uj 22 uj 22 ^ - tr (uj 12 uj 12 ^ + uj 21 uj 21 ^) , 



(B6) 



(B7) 



where in the first equality we have employed the fact that rj is diagonal and that the diagonal elements of uj are 
zero by definition. Each independent variable appears in one and only one of the four terms in the above expression 
so that their contribution to the Berezinian is multiplicative. The contribution from the first term can be read off 
immediately, we obtain 



trg r\ry 



2 kj 

nn% 

j=i P= i 



(B8) 



We now write tru; 11 ^ 11 ^ = J2 P < P ' (| w pp'| 



li 



Each term in the sum contains only independent variables 

which do not appear in any other term so that the contribution of these terms to the Berezinian is multiplicative 
again. Using Eq. (B6) and the anti-hermiticity of du Cl and dv Cl we obtain for the contribution to the Berezinian 



"v. ."I 



fci 



no 

p<p' 



■pi 



s 2 pll ) 2 = A 2 ki (s 2 ) 



(B9) 



In complete analogy, 



,2 \2 
S q'2J 



(BIO) 



q<q' 



Similar arguments are made for the remaining term in Eq. (B7). Since uj 12 and uj 21 couple commuting and anticom- 
muting variables, their contribution to the Berezinian appears in the denominator. Specifically, we obtain 



tr (c/V 2 * + u 2l u> 



21 - 21t )-(flfl(4i+^)) 

\p=lq=l / 



(Bll) 



Collecting terms, we finally obtain the Berezinian (2.26). 

Since Eq. (B3) implies that the metric tensor in the subspace of the radial coordinates is just the unit matrix, the 
radial part of the Laplacian has the form (2.27). 
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APPENDIX C: SEPARABILITY IN SUPERSPACE 



Again, we write A s in a more convenient form, 



We now evaluate the derivatives in Eq. ( |2.29| ) in analogy to the case of ordinary matrices. The calculation is somewhat 
more involved but still reasonably straightforward so that we only mention the intermediate result 

i -3^^CTg|(|; + ^- 4D ™) AW - (C2) 

where D PJ = S p j + s pj (S pj — T p j). Here, S P j = S PJ — S P j and T pj = T PJ — T pj with 

kj ^ kxU) 

Spj = X] „2 _ „2 = X] „2 T „2 (03) 

9=1 ' Pi V g=l Pj qxU) 
q^p 

kj ^ fe xU) ^ 

= X! /-„2 _ „2 \2 ^PJ = X! r« 2 4- q 2 ^2 • ( C4 ) 

,=i y s pj s qj) q =i y s pj + s qx (j) ) 

In the above, we have introduced the convention 

x{j) = {Yi j j Z\ • (06) 

We now show that Y^j=i YlpLi D P j — 0. According to the definition, we have 

D p j — S P j + s p j(S p j — Tpj) — S p j + s p j(S p j — 2S P jS P j + Tpj) ■ (06) 

For each j, the sum over p of the first two terms is zero in analogy to the case of ordinary matrices which was discussed 
in Appendix |A|. We are thus left with the remaining two terms which we denote by —D' p j. Summing over p and j, 
some algebra leads to 



2 kj ki k ? / 



X] D P'J s 2 , 2 

J = l p =l p =l g =l °pl ^ ' <?2 



E 4(4 + s p'l) + 2s pl s p'l ^ s pl( s q2 + s g'2) + 2s g2 S g'2 j ( 
(4i-4i)(4i + 4) h (^-42)(4 + 42) 1 ■ " 



Renaming p <-> p' and 5 <-> in the first and second term, respectively, shows that each of the two sums yields zero 
individually. This completes the proof. 



APPENDIX D: ON BOUNDARY CONTRIBUTIONS IN THE CASE OF HERMITEAN 

SUPERMATRICES 



The purpose of this Appendix is to clarify the role of some boundary contributions to the supersymmetric Itzykson- 
Zuber integral which arise in the case of Hermitean supermatrices. This discussion is not directly related to the main 
content of the present paper. 

After we had computed the result ( |3.30 ) for the function rj(s, t), P.-B. Gossiaux pointed out to us that the structure 
of these contributions in the case of Hermitean supermatrices ought to be very similar. Indeed, this is true. A careful 
reexamination of the considerations following Eq. (B33) in Appendix B of Ref. pleads to additional terms very similar 
to the ones in Eq. (3.30). To clarify this, we calculate the function 



r)(s,t) = 



B k {s) 



det [fj(s p i,is q2 ,t)} 



(Dl) 
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in the case of Hermitean supermatrices for k\ = k 2 — k. Here, the entries of the determinant are given by 

-oo +00 



fj(Spi,is q 2,t) = 



1 

2nt 



dridr 2 
n - ir 2 



exp 



(-^(Gv-^i) 2 + (^2-r 2 ) 2 )^ 



(D2) 



We introduce polar coordinates r\ + ir2 = ftexp(i??) and s v \ + is q 2 — fxexp(iip) and find 



00 2tz 

«*, . •»,». «) - i / * =P (4 (« 2 + / * exp(«) exp (2! 



cos(z9 — -0) 



(D3) 



The angular integration ij4elds the modified Bessel function 7i(/c/x/t). The radial integration can then be performed 
using standard methodsBa Collecting everything, we arrive at 



rj(s,t) 



B k {s) 



det 



1 



Spl - is q2 



1 — exp 



S pl + S q2 

2t 



- p,q— l,...,k 



As in the case of complex matrices, this function satisfies the diffusion equation 

1. , . d . 

-A s T)(S,t) = g^V{s,t) , 



(D4) 



(D5) 



where the radial part A s of the Laplacian is defined in Eq. (B19) of Ref. 0. The main results of Ref. do not change 
due to these additional contributions. 



1 C. Itzykson and J.B. Zuber, J. Math. Phys. 21, 411 (1980). 

2 Harish-Chandra, Am. J. Math. 80, 241 (1958). 

3 A. Pandey and M.L. Mehta, Commun. Math. Phys., 87, 449 (1983); M.L. Mehta and A. Pandey, J. Phys. A 16 2655, L601 
(1983). 

4 M.L. Mehta, Random Matrices, 2nd ed. (Academic Press, New York, 1991). 

5 E. Abdalla, M.C.B. Abdalla, D. Dalmazi, and A. Zadra, ZD-Gravity in Non-Critical Strings, Lec. Not. Phys. m 20 (Springer, 
Berlin, 1994). 

6 S.L. Shatashvili, Commun. Math. Phys. 154, 421 (1993). 

7 T. Guhr, J. Math. Phys. 32, 336 (1991). 

8 F.A. Berezin, Introduction to Superanalysis, MPAM vol. 9, (D. Reidel Publishing, Dordrecht, 1987). 

9 K.B. Efetov, Adv. Phys. 32, 53 (1983). 

10 J.J.M. Verbaarschot, H.A. Weidenmiiller, and M.R. Zirnbauer, Phys. Rep. 129, 367 (1985). 

11 T. Guhr and H.A. Weidenmiiller, Ann. Phys. (NY) 199, 412 (1990). 

12 T. Guhr, Phys. Rev. Lett. 76, 2258 (1996); Ann. Phys. (NY), accepted for publication. 

13 J. Alfaro, R. Medina, and L. Urrutia, J. Math. Phys. 36, 3085 (1995). 

14 T. Guhr, Commun. Math. Phys. 176, 555 (1996). 

M.J. Rothstein, Transactions of the American Mathematical Society 299, 387 (1987). 

16 T. Guhr, Nucl. Phys. A 560, 223 (1993). 

17 J. Ambj0rn, J. Jurkiewicz, and Yu.M. Makeenko, Phys. Lett. B 251, 517 (1990). 
J. Ambj0rn, L. Chekhov, C.F. Kristjansen, and Yu.M. Makeenko, Nucl. Phys. B 404, 127 (1993). 

E. Shuryak and J.J.M. Verbaarschot, Nucl. Phys A 560, 306 (19 93); J.J.M. Verbaarschot, Phys. Rev. Lett. 72, 2531 (1994). 
A.D. Jackson and J.J.M. Verb aarschot, preprin t hep-ph/9509324 , to appear in Phys. Rev. D; A.D. Jackson, M.K. §ener, and 
J.J.M. Verbaarschot, preprint |hep-ph/9602225| . 

T. Wettig, A. Schafer, and H.A. Weidenmiil ler, Phys. Lett. B 36 7, 28 (1996). 
M.A. Stephanov, preprints |hep-lat/960100l| an d |hep-lat/9604003| . 

M.A. No wak, G. Papp, a nd I. Zahed, preprints [hep-ph/9604235| and |hep-ph/9603348|; R. Janik , M.A. Nowak, and I. Zahed, 
preprint hep-ph/960425J ; J. Jurkiewicz, M.A. Nowak, and I. Zahed, preprint |rep-ph/9603308 . 
T. Banks and A. Casher, Nucl. Phys. B 169, 103 (1980). 



15 



L.K. Hua, Harmonic Analysis of Functions of Several Complex Variables in the Classical Domains (American Mathematical 
Society, Providence, RI, 1963). 

G.N. Watson, Theory of Bessel Functions, 2nd ed. (University Press, Cambridge, 1962). 
F. Constantinescu and H.F. de Groote, J. Math. Phys. 30, 981 (1989). 
M.R. Zirnbauer and F. Haldane, Phys. Rev. B 52, 8305 (1995). 



16 



